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ABSTRACT 

We address the effects of bar-driven secular evolution in discs by comparing their 
properties in a sample of nearly 700 unbarred and barred (42 ± 3 per cent of the 
population) massive disc galaxies (M* J? 10 10 M Q ). We make use of accurate structural 
parameters derived from bulge/disc/bar decompositions to show that, as a population, 
barred discs tend to have fainter central surface brightness (A^ ~ 0.25 mag in the i- 
band), and disc scale lengths that are larger by a factor » 1.15 than those of unbarred 
galaxies. The corresponding distributions of fio and h are statistically inconsistent 
at the 5.2 a and 3.8 a levels, respectively. Bars rarely occur in high-surface brightness 
discs, with less than 5 per cent of the barred population having fi < 19.5 magarcsec -2 
- compared to 20 per cent for unbarred galaxies. They tend to reside in moderately 
blue discs, with a bar fraction that peaks at (g — i) « 0.95 mag and mildly declines for 
both bluer and redder colours. These results demonstrate that bars induce noticeable 
evolution in the structural properties of galaxy discs, in qualitative agreement with 
longstanding theoretical expectations. 
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1 INTRODUCTION 

Secular processes are expected to have played a significant 
role in establishing the current properties of massive disc 
galaxies. These are dynamically evolved systems, as evi- 
denced by the presenc e of a population o f discs at z ~ 1 with 



similar scale l ength s l|Lillv et al.lll998l ; ISimard et al.l 



Barden et al. 20051 ) and bar fraction (|jogee et al 



1999 



2004 



Sheth et a.1.1 |200ST ) to what is found in the Local Universe 



This suggests that the quiescent phase of massive discs evo- 
lution started at least 8 Gyr ago, thus leav ing ample time 
for se cular mechanisms to operate (but see iHammer et all 
l2005h . 

Bars perhaps provide the most clear evidence of the 
impact of secular evolution on disc galaxies. Analytical 
and numerical calculations show that they drive a signif- 
icant amount of mass and a n gular momentum redistribu- 
tion in the disc jHohll 1971; Scllwood & Wilk inson! [l99l : 
lAthan assoula 2003; Mar tinez- Valpuesta et af1l2006l ). tunnel- 
ing material towards the galaxy inner regions that can result 
in enhanced gas and stellar den sities - thus possibly fu eling 
nuclear starbursts and AGNs IjShlosman et al.l 1989) and 
leading to the formation of central bulge-like components 
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l|Kormendv fc Kennicuttll2004l ; IDebattista et al.ll2006t ). Fur- 
thermore, the angular momentum exchange between the in- 
ner and outer disc can lead to increased disc scale lengths, 
and the development of ma s s profile breaks at larg e radii 
l|Valenzuela fc Klypinl 120031 : IDebattista et all 120061 ). It is 
precisely the amount of angular momentum exchanged 
within the galaxy what determines the bar strength, and 
this ultimately depends on the mass and velocity distri- 
butions of the material in the disc and sp heroidal (bulge 
plus halo) components l|Athanassoulal [2003) . Once formed, 
these non-axisymmetric perturbations appear to be long- 
lived, surviving buckling instabilities and the growth of a 
significant central mass concentrati on (> 10%, and perhaps 
up to 20% of the or i ginal disc mass; IShen fc Sellwoodll2004l : 



lAthanassoula et al 



gmal i 
l2005l: 



2006). On the 



. IDebattista et al.l 
other hand, Bournaud et al. ( 20051 ) suggest that a substan- 
tial gas component in the disc (~ 7% of the total visible 
mass) might lead to the weakening, or even destruction of 
the bar due to the angular momentum exchange resulting 
from gas inflow. However, this effect is significantly reduced 
when a responsive da rk matter halo model is considered 
l|Berentzen et al.|[2007l ). 

The importance of bars in galaxy evolution studies does 
not only stem from the fact that they can significantly im- 
pact the evolution of a given galaxy, but also because such 
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structures are rather common in disc galaxies. Visual classi- 
fication using photographic pla tes yields a bar fraction of 6 5 
per cent in the RC3 catalogue (|de Vaucouleurs et al.lfl99lh . 
and recent optical studies indicate that approximately half 
of all massive disc galaxies contain bars dMarinova fc Jogeel 
120071 ; iBarazza et atil2008l ; lAguerri et al.ll2009h . The fraction 
of strong bars increases substantially in dust-penetrating 
near-infrared wavelengths, but the total ba r fraction does 
so on l y slightly, as compared to the RC3 (Eskridge et afl 
l200d ; [Menendez-Delmestre et~all 120071 ; iButa et all |2010h . 
It is however now clear that the bar fraction in the Lo- 
cal U niverse is a strong function of the galaxy stellar 
mass dMendez-Abreu et all 120101 ; iNair fc Abraham! l20ld: 
ICameron et al. 2010T ) . and thus consistent results can only 
be obtained when comparing samples well matched in M*. 
Even though wavelength coverage, detection technique, and 
-possibly to a larger extent- sample selection, all play a role 
in the detection of bars in disc galaxies, most studies point 
to fractions ~ 50 per cent in L > L* galaxies. 

The picture that emerges from theoretical work is cor- 
roborated, at least qualitatively, by a number of observa- 
tional results. The bar-driven redistribution of angular mo- 
mentum affects the interstellar medium in galaxies, resulting 
in flatter chemical abundance gradients in barred galaxies 
|Martin fc Rovlll994l ; rZaritskv et al.lll994l ), as well as higher 
central concentrations of molecular gas (|Sakamoto et alj 
1 19991 ; ISheth" et al. I l2005l ). More recently, evidence has been 
found that the gas brought t o the centre b y bar s is effi- 
ciently transformed into stars. [Ellison et all 1 201ll) present 
indication that the current star formati on rate at the cen- 
tre is higher in massive barred galaxies. ICoelho fc GadottH 
(2011) show that the distribution of mean stellar ages in 
bulges of massive barred galaxies shows a peak at low ages 
that is absent for their unbar red counterparts (see also 
IPerez fc Sanchez-Blazauezll201ll ). 

Yet investigation of the effects of bar-driven secular evo- 
lution on the structural properties of discs is currently lack- 
ing in the literature. In this Letter we fill this gap by com- 
paring the disc properties in sample of nearly 700 barred 
and unbarred galaxies. 



2 GALAXY SAMPLE 



The s ample used here is the one presented in iGadottil 
(|2009h n. which contains all galaxie s in the SDSS-DR2 wit h 
stellar masses M* > 10 10 M ffrom lKauffmann et al. 2003^ 
at redshift 0.02 ^ z ^ 0.07, and with axial ratio b/a ^ 0.9-13 
These criteria provide a sample which is both representative 
and suitable for 2D bulge/disc/bar decomposition. The red- 
shift range allows for enough spatial resolution, while select- 
ing face-on galaxies minimizes dust and projection effects, 
and eases the identification of bars. The reader is referred 
to that paper for a detailed discussion of selection effects. 

Through multi-band (gri) 2D decomposition, IGadottil 

(2009) provides accurate structural parameters for the three 
components (when necessary) - including the central surface 
brightness, scale length, integrated colour and stellar mass of 



1 See http://www.sc.eso.org/~dgadotti/buddaonsdss.html 

2 Throughout this paper we assume a flat cosmology with Qjvi 
0.3, Q A = 0.7 and H Q = 75 kms" 1 Mpc" 1 . 



the disc component, as well as the galaxy bulge-to-total ra- 
tio. Despite of the inhe r ent co mplexity of bulge/disc/bar de- 
compositions, Gadotti (2009, Appendix A) shows that disc 
structural parameters are particularly stable, and therefore 
both h and no can be robustly derived even in the presence 
of a bar. Moreover, it is important to recall that failing to 
account for the contribution of the bulge results in structural 
parameters corresponding to a maximum disc, and therefore 
our multi-component fits provide the least biased results. In 
this Letter we focus on structural parameters derived from 
fits in the i-band, as it is the least affected by dust and re- 
cent star formation. Disc stellar masses were obtained from 
disc luminosities and mass-to-light ratios in the i-band, the 
latter derived from the integrated (g — i) disc co lour and the 
relation determined bv lKauffmann et al.l (|2007l ). In order to 
select a clean sample of disc galaxies, only systems having 
B/T < 0.8 have been considered. 

To verify whether a galaxy is barred, typical bar sig- 
natures were searched for through the inspection of each 
galaxy image, isophotal contours and a pixel-by-pixel ra- 
dial intensity profile. It should be noted that due to the 
limited spatial resolution of SDSS images we miss most 
bars with semi-major axis shorter than Lbar ~ 2kpc, 
whi ch are mainly found in very lat e-type spirals (later than 
Sc; lElmegreen fc Elmegreenl 1 19851 ) an d are usually not de- 
tected in most re c ent stud i es (e .g., IBarazza et al.l 120081 ; 
lAguerri et al.l200lil ). IGadottil (|201lf ) presents a detailed anal- 
ysis of the structural properties of the bars in this sample. 

Our final sample consists of 291 barred and 393 un- 
barred disc galaxies, corresponding to a bar fraction of 42 ± 
3 per cent (binomial 90% uncert ainty) - in excellent agree- 
ment with previous studies (e.g.. lAguerri et "ah 1 120091 ). Both 
subsamples are perfectly matched in disc stellar mass, with 
2.5, 50 and 97.5 per cent quantiles of [0.5,1.7,5.5] x 10 10 M Q 
and [0.4,1.6,6.5] x 10 10 Mq for barred and unbarred galaxies, 
respectively. 



3 PROPERTIES OF DISCS IN BARRED AND 
UNBARRED GALAXIES 

Figure[T]reveals the striking differences of the structural scal- 
ing relations of galaxy discs in barred and unbarred galaxies. 
Panels la and 16 show, respectively, the i-band disc central 
surface brightness (/xo) and scale length (h) as a function of 
disc stellar mass0, while panel lc depicts the bivariate distri- 
bution of the two former quantities. Isocontours in all panels 
enclose 25, 75 and 95 per cent of each population. There are 
several remarkable features worth discussing. First, panel la 
shows that the distribution of disc central surface brightness 
is significantly different for barred and unbarred galaxies, in 



3 By using the current disc stellar mass instead of the galaxy 
total stellar mass we eliminate all the possible uncertainties and 
biases associated to varying bulge-to-disc ratios at fixed M 4 . We 
also constructed this plot using the bar-plus-disc mass in the case 
of barred galaxies, to represent the original disc mass before the 
instability developed. This does not modify our results, simply 
because the bar mass fractions are always small (median bar-to- 
total ratio ~0.1; IGadottil l201lf) . We however prefer to use the 
current disc mass as derived from the exponential fits because it 
is the quantity directly related to h and fig. 
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Figure 1. Structural scaling relations of discs in barred and 
unbarred galaxies. In all panels isocontours enclose 25, 75 and 95 
percent of each population, a) Disc central surface brightness, /iq, 
as a function of disc stellar mass, b) The disc scale length vs stellar 
mass relation, c) Joint fio-h relation for both subsamples. All 
parameters are derived from fits to the i-band images. It is clear 
that discs in barred galaxies are characterised by having fainter 
central surface brightness and larger scale lengths. Note the lack 
of bars in galaxies having compact, high surface brightness discs. 
Shaded regions indicate discs where we likely miss more than half 
of the bars because of their small length. 



the sense that at a given disc stellar mass barred galaxies 
tend to have fainter fio values - note how all isocontours 
extend towards much brighter /xo for unbarred galaxies than 
for barred ones. Notably, the reduction in surface brightness 
is accompanied by an increase of disc scale length (panel 
16), resulting in markedly dissimilar distributions in the /io- 
h plane (panel lc). While unbarred galaxies populate the 
high-surface brightness and small-scale lengths parameter 
space, barred galaxies are essentially absent from this re- 
gion - only 5 per cent of all barred galaxies have no < 19.5 
magarcsec -2 , while this fraction increases up to 20 per cent 
for unbarred discs. One potential explanation for the lack of 
bars in these discs is that they have semi-m ajor axis l ength s 
Lbar < 2 kpc and thus remain undetected. iGadottil (|201 lh 
shows that there is a tight relation between Lf, ar and h, 
with a median ratio Li, ar /h — 1.5 that appears to be inde- 
pendent of disc scale length within this sample. This implies 
that we likely miss more than half of the bars in discs having 
h < 1.33 kpc. The shaded areas in Fig.[T]highlight this region 
of incompleteness, and show that the paucity of bars in high 
surface brightness discs is real. First, we still detect bars in 
these small discs, but they only occur in the lowest surface 
brightness systems. Second, panel lc shows that bars are still 
missing in larger discs having fio < 19.5 magarcsec -2 , where 
bar non-detection is not an issue. The remarkable result is 
the paucity of long bars in these systems, and the evidence 
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Figure 2. Normalised distributions of disc central surface bright- 
ness (a) and disc scale lengths (b) for barred and unbarred galax- 
ies. According to a Kolmogorov-Smirnov test, the corresponding 
distributions of fio and h are statistically inconsistent at the 5.2 a 
and 3.8 a levels, respectively. As a population, discs in barred 
galaxies tend to have £«0.25 mag fainter central surface bright- 
ness and ?s 15 per cent larger disc scale lengths. 



that barred discs extend towards the large-fr, faint-/io region 
of the plot. 

These differences can be more clearly appreciated in 
Fig. [2] where we show the normalised distributions of fio and 
h for our two populations (panels 2a and 2b, respectively). 
The central surface brightness distribution of discs in un- 
barred galaxies peaks at a moderately more luminous value, 
and features a distinct extended tail towards brighter fio- 
A Kolmogorov-Smirnov test rules out the null hypothesis 
that the two distributions are drawn from the same parent 
population at the 5.2 a level. The distributions of disc scale 
lengths for both populations are also inconsistent at the 3.8 a 
level according to the KS test. The two distributions would 
however be statistically indistinguishable if barred galaxies 
had 0.25 mag brighter central surface brightness and 15 per 
cent smaller scale lengths. It is important to note that, de- 
spite the high statistical significance of the differences here 
reported, there is considerable scatter at a given disc stellar 
mass: cr M0 ~ 0.5 mag and an ~ 1.10 kpc. 

Aside from the structural scaling relations, we investi- 
gate the stellar population properties of discs in our sample 
of barred and unbarred galaxies. Fig.[3ji presents the rela- 
tion between the disc stellar mass and its integrated (g — i) 
colour, while panel [3j) shows the corresponding normalised 
colour distributions. The disc colour distribution appears to 
be marginally bimodal for both populations^ but, remark- 
ably, their peaks differ from (g — i) ~ 1.25 for unbarred 
galaxies to (g — i) ~ 0.95 for barred ones. The dark gray 
curve in Fig.[3J) shows the bar fraction as a function of disc 
colour (only for bins having more than 5 galaxies), while the 
light gray curves correspond to the binomial 90% confidence 
intervals. Bars tend to reside in moderately blue discs, with 
a fraction that peaks at ~ 55% at the previously mentioned 
value and mildly declines for both bluer and redder discs - 
in qualitative agreement w ith previous work l|Aguerri et al.l 
120091 : iMasters et~aH l201lT ). The fraction of bars in galax- 



A one-sided Kolmogorov-Smirnov test rules out the null hy- 
pothesis of a normal distribution at the 90% (1.6 a) and 99.8% 
(3.1 cr) levels for barred and unbarred galaxies, respectively. 
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Figure 3. a) The disc (g — i) colour vs stellar mass relation for 
barred and unbarred galaxies, b) Normalised distributions of disc 
colour for both subpopulations. The dark gray curve shows the 
strong dependence of the bar fraction with disc colour, while light 
gray curves correspond to the binomial 90% confidence intervals. 
Bars tend to reside in moderately blue discs, with a bar fraction 
that peaks at 55% at a disc colour of [g — i) Ri 0.95 mag, and 
declines for both redder and bluer colours. 



ies with very blue col ours is probably a lower limit (cf. 
iNair fc Abrah^ml I2010T I. as we recall that these late-type 
systems can contain small bars that are undetected in our 
images. 



4 DISCUSSION AND CONCLUSIONS 

We have shown that barred and unbarred M* > 10 10 M© 
galaxies are characterised by distinct distributions of disc 
structural parameters. As a population, barred discs tend to 
have fainter central surface brightness (A/io ~ 0.25 mag in 
the i-band) and disc scale lengths that are larger by a factor 
~ 1.15 than those of unbarred galaxies. This value is close 
to, but slightly lower than, that predicted by the numerical 
simulations of lValenzuela fc Klypinll|2003l ), where discs were 
found to increase their scale lengths by factors 1.2—1.5 due 
to the transfer of angular momentum that accompanies the 
formation of a bar. Nevertheless, the moderate increment 
of scale lengths revealed by the data has to be thought of 
as a lower limit to the actual effects of bar-driven secular 
evolution. The reason is twofold. 

First, it is important to recall that bars are not the 
only drivers of disc secular evolution. In fact, any type 
of non-axisymmetric perturbation -including spiral arms, 
oval distortions and triaxial dark matter haloes- can also 
modify the propertie s of discs (e.g- jKormendv fc Kennicuttl 
|2004| ; [Sellwoodll2010l ). but they generally operate on longer 
timescales. It is therefore possible that secular evolution has 
also occurred, to some degree, in our sample of unbarred 
galaxies. 

Second, our disc structural parameters come from fit- 
ting one single exponential to the surface brightness profile, 
with no distin ction between different disc Freeman types 
(Freeman 1970). This can slightly bias our recovered param- 
eters, resulting in marginally brighter fio and smaller h in the 
case of Type II profiles with outer truncations - which occur 
far more often than sin gle exponential, Type I profiles (e.g., 
iPohlen fc Truiilloll2006l ). This is indeed supported by the re- 



cent analysis of disc properties in S 4 G galaxies l|Sheth et alj 
2010). Munoz-Mateos et al. (in preparation) find that, when 
allowing for outer profile breaks, disc scale lengths in barred 
galaxies are, on average, a factor ~ 1.8 larger than those 
without bars - but, as in our case, the scatter at fixed 3.6 /im 
magnitude is large. 

In any case, our fits provide the optimal average expo- 
nential profile between small and large radii, and thus allow 
for a direct comparison with the numerical simulations by 
iDebattista et alT(|2006l ). Not surprisingly, they find that the 
amount of evolution in their simulated discs depends criti- 
cally on the initial disc kinematics. The density profiles of 
highly unstable discs (low Toomre's Q) evolve dramatically 
compared to more stable initial configurations, resulting in 
final scale length differences of factors > 2 even for models 
with nearly identical initial angular momentum. The gen- 
eral situation is of course more complex than this, and the 
specific density profile evolution in their simulations is deter- 
mined by the phase-space distributions of the stellar disc and 
the dark matter hal o - resulting in scale lengths c hanging by 
factors 1.0-2.4. As IDebattista et al.l (fcOO-jfeoOfj ) point out, 
this has the important consequence that direct estimates of 
dark matter halo spin parameters from measured disc scale 
lengths can be rather uncertain. 

Our results support their claims, but suggest that the 
amount of scale length evolution due to bar formation is 
only moderate. A simple back-of-the-envelope argument ap- 
pears to be consistent with this idea. If secular mechanisms 
were to increase disc scale lengths by a considerable amount 
-say, factors larger than two-, we would then expect to see 
clear signs of evolution in the mass-size relation of discs be- 
tween redshifts < z < 1. Yet observations support mild 
to no evolution at all (|Lillv et al.lll998l ; ISimard et al.lll999l ; 
iBarden et ai1l2005h and this, in turn, is consist ent with real 
discs being reasonably stable (Q > 1.5; e.g., iKregel et al.l 
2005). If this is the case, and considering all the assump- 
tions involved, it is most likely that the uncertainties de- 
rived from mapping halo spin to disc scale lengths are not 
dominated by secular evolution effects - but they certainly 
are a contributing factor. 

It is not yet clear if there exists one single condi- 
tion d^be^mfriing whether a galaxy will be bar-stable or 
not. lAthanassoulal (|2008T ) shows that simple criteria that do 
not fully capture the complexity underlying bar formation 
[e.g., the Efst athiou et alj l| 19821 ') criterion] generically fail 
to correctly predict disc stability. Instead, numerous fac- 
tors come into scene in order to stabilise a disc against 
non-axisymmctric perturbations. Thus, galaxies with very 
weak or no bars must have either a kinematically hot 
disc and/or a significant central mass concentration and/or 
a very low relative disc mass and/or be embedded in a 
quite unrespon sive dark matter halo ( Athana ssoulal 120031 ; 
ISellwoodl |2010| ) . From the disc component point of view, 
a high stellar velocity dispersion can provide significant 
stabilisation an d prevent the formation of a b ar for over 
a Hubble time llAthanassoula fc Sellwoodlll986l ). Moreover, 
ISellwood fc Evans! ( 200lh show that a steeply rising inner 
rotation curve is sufficient to bar-stabilise a disc, regardless 
of the dark matter content. In this context, the realisation 
in Fig.[]J that compact, high-surface brightness discs do not 
host bars -or if they do, they are rather small- is intrigu- 
ing, and it is tempting to think of these systems as having 
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a high rate of shear at the centre capable of (almost) fully 
stabilising the disc. Kinematical studies of these galaxies are 
highly desirable in order to test this hypothesis and clarify 
their lack of bars. 

Finally, basic stability criteria and numerical simula- 
tions suggest that bars are more easily triggered in gas-rich 
(i.e., cold), star-forming discs. This provides a natural expla- 
nation for the preferential occurrence of bars in moderately 
blue discs, but the exact shape of the bar fraction distri- 
bution deserves further investigation. Even though such an 
analysis is beyond the scope of this Letter, one can naively 
think that gas consumption through star formation leads 
to redder and kinematically hotter (more bar-stable) discs. 
In this context, the declining fraction towards bluer colours 
has probably to be understood in terms of a selection effect, 
such that these galaxies contain smaller bars that remain 
undetected in our images. 

We have shown that bars drive noticeable structural 
evolution of galaxy discs, confirming longstanding analyt- 
ical and numerical predictions. Detailed structural decom- 
position of galaxies provides one of the most powerful di- 
agnostics for galaxy evolution studies, allowing for direct, 
quantitative comparisons with theory and simulations. 
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